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This talk is about solving cosmological equations analytically without approximations, and dis- 
covering new phenomena that could not be noticed with approximate solutions. We found all the 
solutions of the Friedmann equations for a specific model, including all the zero-size-bounce solutions 
that do not violate the null energy condition, as well as all the finite-size-bounce solutions, and then 
discovered model independent phenomena. Among them is the notion of geodesic completeness for 
the geometry of the universe. From this we learned a few new general lessons for cosmology. Among 
them is that anisotropy provides a model independent attractor mechanism to some specific initial 
values for cosmological fields, and that there is a period of antigravity in the history of the universe. 
The results are obtained only at the classical gravity level. Effects of quantum gravity or string 
theory are unknown, they are not even formulated, so there are new theoretical challenges. 



I. INTRODUCTION 



This work started with the discovery of new techniques for solving the Friedmann equations analytically [1| , 
and finding the complete set of solutions for the gravitational system that includes a scalar field a (x^) with 
a special, but still typical, potential V (a) in a cosmological model [H0]- The special techniques were direct 
outcomes from 2T-physics 0] and 2T-gravity This led to an emergent formulation of gravity in lT-physics 

j— i' with an additional local conformal (Weyl) symmetry. The cosmological model that emerged in this approach 
had a strong overlap with the cosmological colliding branes scenario 6] 7} developed from worldbrane notions 
Q in M-theory Q. This connection led to a fruitful collaboration between our group, that includes C.H. Chen, 
5h , and the group that includes Paul Steinhardt and Neil Turok. Our collaborative work has partially appeared in 
, [To| . Here I will summarize results that include some that are about to be published in fill ] . 

Our approach begins with the standard action typically used in cosmological models that describe a scalar 
CN| field a minimally coupled to gravity 

<N i s= f d i xV - f (ff) -kf^a- V(a)\ 

/ ^radiation + matter I 

00 ■ } 

Although not necessary, the scalar field a could be identified with the dilaton in string theory, written in 
ON ' the Einstein frame. For the purpose of studying cosmology, the relevant degrees of freedom are homogeneous 
fields that are only time dependent. The space dependence is ignored in the leading approximation. It can be 
argued that the space gradients are responsible for non-leading effects as compared to the dominant effects 
. . , parameterized by only time dependent fields that we study here (including anisotropy). 

Our cosmological analysis of this theory involves geometries with or without spacial curvature, as well as 
anisotropy, and includes radiation and other relativistic matter approximated as relativistic dust described by 
an energy-momentum tensor consistent with conformal symmetry. The geometries we study have the form 



ds 2 = -dt 2 + a 2 (t) dsl = a 2 (r) (-dr 2 + dsf) , (2) 



where a 2 (r) is the scale factor that tracks the size of the universe. We will use conformal time r which is related 

4 



to cosmic time t, by dt = a (r) dr, as in Eq.([2]). The 3-dimensional part ds 2 of the geometries we study include 



the isotropic FRW universe with or without curvature K 



(dsl) FRW = T~T7~a + r 2 (dO 2 + sin 2 Od^ 2 ) ; K = * (3) 



n 







3 J FRW I _ Xr 2 

• the anisotropic Kasner metric with zero curvature 

( d 4) Kasner = e- 2 VV~^ {dz f + e V57W ( dx ) 2 + e -V5«» {dy) 2 ) , (4) 



where oti (r) and cti (r) are two dynamical degrees of freedom in the metric along with a (r) , and 
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• the anisotropic Bianchi IX metric with curvature K 

{dsl) IX = e~ 2 VV~^ (da z f + e VV~^ h^™* (da x ) 2 + e"^" 2 (da y ) 2 ) , (5) 

where d<Ji, which are given in [liij . depend on space coordinates (and K) but not on r. The Bianchi IX 
metric reduces to the FRW metric with curvature K in the isotropic limit a* — > 0, and it reduces to the 
anisotropic Kasner metric when the curvature vanishes K — > 0. 

The relevant degrees of freedom are a (r) , (r) , a\ (r) , «2 (t) . Here we have attached the label £7 on the 
scale factor a e (t) to emphasize that it is defined in the Einstein frame. We do this because we will also 
define other frames to solve the equations. The equations of motion (namely the Friedmann-type equations) 
are obtained from the action above. However, it is convenient to provide an effective worldline mechanics-type 
action in the Einstein frame S^ s that reproduces the same equations 

S ^-J dT { -e[a E V{a)+p -^alv{a u a,)} )' (6) 

Here the constant po represents the energy density of radiation and relativistic matter. It comes from approxi- 
mating these degrees of freedom as "dust" modeled as a perfect fluid, with an energy-momentum tensor of the 
form Tfu; = (p + p)u ll u u +pg[_iv, that is taken to be traceless (p = 3p) and behave as p(r) = pa/a% (r) to be 
consistent with the conformal symmetry for massless particles (including radiation and matter). The potential 
energy term for the anisotropy involving v {ct\, 02) 

v («!, a 2 ) = ^ (-e-VaTW - 4e 2 V^^i sin h 2 ( V2^a 2 ) + 4e~VV^ Ql cosh (V2 K a 2 ) ) (7) 

is obtained by computing R (g) for the Bianchi IX metric [13( . The additional degree of freedom e (r) plays 
the role of the metric on the worldline and ensures r-reparameterization invariance. This effective action 
can be derived directly from the 4-dimensional action by inserting the metric including the lapse function 
ds 2 = a? E (r) (— dr 2 e 2 (r) + ds 2 ^j , where the lapse function e (r) is the einbein on the worldline. The equation 
of motion for e (r) imposes the zero energy constraint which is equivalent to the Goo = Too Einstein equation. 
Eventually we work in the gauge e (r) = 1 which corresponds to the conformal time in Eq. ([2]) . Similarly, it can 
be checked that this effective action reproduces all Friedmann-type equations that are derived from the original 
action in Eq.JT]). 

For our model we chose the following potential energy 

V{a) = (A)' (csinh 4 U£*\ +6cosh 4 (j£ ajj , (8) 

After insuring stability by requiring (b + c) > 0, the plot of this potential, which looks like a well or double 
well for various values and signs of 6, c, shows that it has features that are similar to other potentials used 
in cosmological applications. The reason for choosing this specific potential is that we can use some special 
tricks to solve the Friedmann equations exactly in the FRW universe and find all of the solutions analytically 
[H 03 03 • We can similarly solve the equations analytically for a few other potentials, as presented in a future 
publication. 

II. BCST TRANSFORMATION AND WEYL SYMMETRY 

The tricks for solving the equations are actually deeply related to the gauge symmetries and shadow phenom- 
ena in 2T-physics [5|[l|. But, in the current setting with the potential V (a) in jSJ), it is possible to introduce 
them directly as the Bars-Chen-Steinhardt-Turok (BCST) transformation of the fields (a%,oj to a new basis 
(4>j, s 7 ) as follows 
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Then the action S^L in Eq.© is rewritten in the new basis as follows 



Q7 . 
D cff ' 



dr 



is 2 

2*7 



an 



[0 7 /(s 7 /0 7 ) +Po ~ f (07 - s 2 J )v(a 1 ,a 2 )] 



(10) 



cs 7 when V (a) is given as in Eq.([8]). 



where tf^f (s 7 /</> 7 ) = &</> 7 , . . 

In the absence of anisotropy, cu —> 0, u (ai, «2) — > 1, this form of the action, with the special potential, shows 
that the fields (</> 7 ,s 7 ) are decoupled from each other, except for the zero energy condition that results from 
the e-equation of motion. The decoupled second order field equations for <f> 1 (r) and s 7 (r) are easily solved 
because they are analogous to a particle in a potential. The solutions are given in terms of Jacobi elliptic 
functions. The zero energy condition merely relates the integration variables E^, E s (energy levels for 7 ,s 7 ) 
in the form E^ = E + po and E s — E, so that there is only one integration variable, E, and two initial values 
7 (to) , s 7 (to) . Due to time translation symmetry one of the initial values has no physical significance; therefore 
E,4>{jo) are the only two integration parameters. Together with the 4 parameters of the model (&, c, K, p$) , 
these 6 parameters determine the properties of all the solutions in the absence of anisotropy. 

This shows that the BCST transformation makes the Friedmann equations solvable analytically and leads to 
the complete set of solutions given in [l[ 1(| U| as outlined below. 




Fig.l- Gravity (f> 2 > s 2 in left/right quadrants. Antigravity <fi 2 < s 2 in top/bottom quadrants. 

Eq.Q is only half of the BCST transformation between the (a%, a) and (0 7 , s 7 ) bases. The form in Eq.© is 
initially defined only for the region of the (<^ 7 , s 7 ) plane that satisfies (4> 2 — s 2 ) > 0, which corresponds to the 
left and right quadrants as depicted in Fig.l. The cosmological singularity at a 2 E = corresponds to the 45° 
lines, which may be considered a "lightcone" in the (</> 7 ,s 7 ) space. The study of geodesic completeness at the 
a% = singularity amounts to the question of how geodesies cross the lightcone in (</> 7 , s 7 ) field space. Through 
our explicit analytic solutions [l| [l(| [ll| we learned that generic trajectories in the absence of anisotropy cross 
this lightcone and continue to develop in the entire (0 7 , s 7 ) plane. Therefore, geodesic completeness leads us to 
expand the domain of the Einstein frame beyond the (</> 7 — st) > region by including the full (0 7 , s 7 ) plane. 
This extension of the field domain, including the top and bottom quadrants, reminds us of a similar extension 
of spacetime (rather than field space) via the Kruskal-Szekeres coordinates of the Schwarzschild blackhole. 

To complete the BCST transformation, we also need the inverse transformation which involves the four 
quadrants in the (</> 7 , s 7 ) space. This is given by 

f fVW\^M , « , > o f |^s inh , » , > o 

1 &yM<to>>{%). i! *«> \ ^v1^i<:™h(^), if z < 

Actually the cj) 2 > s 2 and 4> 2 < s 2 field domains are separated from each other by the spacetime singularity 
at a 2 E = 0; so cf> 2 > s 2 is satisfied in a patch of spacetime, while cf> 2 < s 2 is satisfied in a different patch of 
spacetime. These patches communicate with each other through the cosmological singularity, or equivalently 
through the whole lightcone in the (0 7 , s 7 ) plane. 

Allowing (</> 7 , s 7 ) to be extended to the full (</> 7 , s 7 ) plane, including the top and bottom quadrants (as guided 
by the solutions), defines a continuation of the Einstein frame to a larger domain of field space. We found that 
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the full (A 7 ,s 7 ) provides a geodesically complete geometry, while the domain A 2 > s 2 by itself is geodesically 
incomplete. We find that in the extended domain the physics described by S^g includes patches of spacetime in 
which there is gravity when A 2 > s 2 , or antigravity when A 2 < s 2 . The actions S^ s , are equivalent to each 
other via the BCST transformation only in the sector in which yA 2 — sz) is positive. In the antigravity patches 
the inverse BCST transformation produces again an Einstein frame, and an action similar to S^ s , but with a 
negative Newton constant. 

To better understand the meaning of the different patches we need the complete theory, not only the effective 
theory S7a for the cosmological degrees of freedom. This is given from first principles by a slight extension of 
the Einstein action by including only gauge degrees of freedom associated with a local scaling symmetry, called 
Weyl symmetry. The local scaling transformation, given by 

(<f>(x),s(x))-> (A(x),s(x))e x ^ and 9fat {x) -> 9fiu (x) e~ 2X ^ , (12) 

does not allow the usual Einstein-Hilbert term, but allows conformally coupled scaiars. The following action 
contains two conformally coupled scaiars A, s, interacting with the curvature term with the coefficient dictated 
by the gauge symmetry 

S = J d^xj—g f^d^A - Ig^dpsfa + ^ (A 2 - s 2 ) R (g) - A\f f^j ) . (13) 

Here f(z) is an arbitrary function of the scale invariant ratio z = 4. Fermionic and gauge fields, as well as more 
conformally coupled scaiars, can be added, as part of a complete model, without breaking the conformal gauge 
symmetry [5J. Here (A 2 — s 2 ) /12 appears as a gravitational parameter that replaces the Newton constant 1/2k 2 . 
To have the possibility for it to be positive, we must have one scalar, namely A, to have the wrong sign kinetic 
term. If A, like s, has the correct sign, then we would end up with a purely negative gravitational parameter. 
Hence the relative minus sign in (</> 2 — s 2 ) is required. The wrong sign kinetic term makes A potentially a ghost. 
However, the local gauge symmetry compensates for the ghost, thus insuring unitarity of the theory [19[. 

The usual Einstein gravity in Eq.([T]) can be obtained by choosing a gauge for the local scale symmetry. 
In particular, in the Einstein gauge we denote the fields by A E (x) ,s e {x) , g E (x) with a subscript E, and we 
choose ttj (cj) E (x) — s 2 E (x)) = for all spacetime thus generating the gravitational constant. Then writing 

A E (x) = ±^ cosh (^M) , SE (x) = ±^ sinh (^M) , (14 ) 

we find that the conformally gauge invariant action (|13|) yields the familiar action in Eq.([T]) with an Einstein- 
Hilbert term -^iR (se) and a minimally coupled scalar field a, and any potential energy V (a) . 

Another convenient gauge choice is the 7-gauge in which we denote the fields by <f> 1 (x) , s 1 (x) , g^ u (x) with 
a subscript 7. In this gauge we choose the conformal factor of the metric a 2 (x^) = 1 for all spacetime {x^} . 
In particular, in the FRW, Kasner or Bianchi IX spacetimes of Eqs.(|2][5|) taken in the gamma gauge a 2 (r) = 1, 
we obtain simple expressions for the curvature, such as R [g^ RW ) = 6K, etc.. Inserting these expressions in 
the full action in Eq. (fl~3|) we derive the gauge fixed action 5j ff for cosmological applications, precisely as given 
in F.q.lfTOl). 

It is useful to point our that the following quantity is gauge invariant in all spacetime 

z{x^ = ^{-g{x)^{6 2 {x)-s 2 (x)). 

In the Einstein gauge in the gravity sector given in Eq.fTTJJ it is equal to z (x^) = {—g (x)) 1 = a% (x), while 
in the Einstein gauge in the antigravity sector (interchanging cosh and sinh in Eq.([T4])) it is equal to z (x^) = 

— (— g (x)) J = —a E (x) . Furthermore, in the 7-gauge it is equal to z (x M ) = ^ (</> 2 — s 2 ). Equating these 
quantities to each other by gauge invariance immediately gives part of the BCST transformation in Eq.®. 
This shows that the BCST transformation in Eqs. (|9lll[) is actually the gauge transformation from the Einstein 
gauge (both gravity and antigravity regions) to the 7-gauge and vice- versa. 

Now it is clear that the 7-gauge of the action (flU)) is valid in the entire (</> 7 (x) , s 7 (x)) plane in all spacetime, 
for all values of the gauge invariant z (x). But the usual Einstein theory as given in Eq.([T]) is possible only in 
the left and right quadrants of Fig.l where (<^> 2 — s 2 ) > 0, or equivalently when the gauge invariant z (x) is 
positive. To give an Einstein frame description of the top and bottom quadrants of Fig.l, where (<fi 2 — s 2 ) < 
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and the gauge invariant z (x) is negative, we need an Einstein gauge just like Eq. (fT4"|) by interchanging cosh 
with sinh . The result of gauge fixing in the region (0 2 — s 2 ) < is an Einstein action just like (JT]), but with the 
gravitational constant 1/2k 2 replaced by — 1/2k 2 , and the wrong sign kinetic term for the sigma field. Evidently, 
in this system there is antigravity rather than gravity since the Newton constant has the opposite sign. 

Comparing the gravity/antigravity regions further in their respective Einstein frames, note that the potential 
energy V (a) does not change sign, but is replaced by a new one V (er) which is related to the old one in Eq.® 
by interchanging cosh with sinh . Any additional terms in a complete action do not change sign either, but any 
terms that initially contain <j), s would look different in the gravity/antigravity sectors by the interchange of 
cosh with sinh. The 7-gauge interpolates smoothly between the gravity and antigravity regions and this is why, 
among other things, it is helpful with the geodesic completeness of the geometry. 

After discussing the classical solutions of the cosmological equations, we will return to the antigravity region 
to discuss possible instabilities due to the sign changes in the full action (|13p. and what they may mean. 

III. PROPERTIES OF THE SOLUTIONS 

First, recall the geodesic incompleteness of the geometry in the usual Einstein frame, which is an old problem 
that has been set aside. This is typical for any FRW cosmological solution with a big bang singularity, because 
the geometry is incomplete at the singularity as a% (r) —> 0. A massive particle in the flat FRW universe 
satisfies the geodesic equation dx/dr = pj ^Jp 2 + m 2 a 2 E (r) where p is the conserved momentum of the particle. 
From this it is straightforward to compute that it takes a finite amount of conformal time, cosmic time, or 
proper time, to reach the big bang singularity from anywhere in the universe. Therefore the particle trajectory 
is artificially stopped at the singularity in a finite amount of proper time. This is geodesic incompleteness. It 
begs the question: what was there just before the big bang? What would we find if we could allow our watch to 
run beyond the singularity where it was artificially stopped? A geodesically complete geometry would provide 
at least a partial answer to this question. 

Of course, modifications of the geometry or even a completely different description may be expected from 
quantum gravity or string theory. Unfortunately, the question remains even more obscure in these formalisms 
because of the lack of proper understanding of how they should be applied to the cosmological setting. It is 
therefore not unreasonable to at least, at first, try to answer such a question in the setting of classical gravity. 
Having the complete set of solutions analytically allowed us to answer such questions, and understand how to 
complete the geometry in the classical physics setting, thus leading to some surprising behavior, as follows. 

The complete set of solutions of the cosmological equations have been analyzed in detail in [l[ [H[ ■ These 
solutions are obtained with no restrictions on either the parameters of the model or the initial conditions on the 
fields. Before the important role of anisotropy is taken into account near the singularity at a\ = 0, there are 
6 unrestricted parameters in our solutions as indicated following Eq. (|I0[) . All 6 parameters are available to try 
to fit cosmological data far away from the singularity. 

In the absence of anisotropy, the generic solution for (0 7 (r) , s 7 (r)) behaves in various detailed ways in 
various regions of the 6 parameter space. This is given in the Appendix of [II] where we identified 25 different 
regions of the 6 parameter space in which the analytic expression is different for each case separately. The 
trajectory of the generic solution in the absence of anisotropy can be plotted parameterically in the (0 7 ,s 7 ) 
plane of Fig.l. This shows that the generic trajectory crosses the lightcone in Fig.l at any point (as determined 
by the 6 parameters) so that <fi 2 (r) — s 2 (r) keeps changing sign as the universe moves back and forth from the 
gravity patch to the antigravity patch. Of course, to do so, the scale factor o? E (r) vanishes at each crossing, so 
that the gravity/antigravity patches are connected to each other only through the cosmological singularity. 

The behavior of the generic solution near the singularity changes dramatically in the presence of anisotropy. 
As will be outlined below, an important effect of anisotropy is to focus the trajectory of the generic solution to 
pass through the origin of the (</> 7 ,s 7 ) plane, such that the generic trajectory can cross the lightcone in Fig.l 
only at the origin where (</> 7 ,s 7 ) vanish simultaneously. 

A. Zero-size and Finite-size Bounces in the Gravity Patch 

By restricting the parameter space or initial conditions, we find that there is a subset of special solutions that 
are geodesically complete in the Einstein frame purely in the gravity patch. These are the zero-size and finite- 
size bounces shown in Figs. (2,3,4,5). As seen in Figs. (4,5), their trajectories never reach into the antigravity 
sector and they never cross the lightcone in Fig.l except at the origin. 
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Fig. 2 - Zero-size bounces in a cyclic universe. Fig. 3 - <\> (r) , s (r) for b < 0, with synchronized 
cle (t) and a (r) for b < is depicted. For b > initial values and relatively quantized periods, 
turnaround is at cle — > oo. For b > turnaround is at (f> — » oo. 



The zero-size-bounce (see Figs. (2, 3)) describes a cyclic universe that contracts to zero size and then bounces 
back smoothly from zero size. It expands up to either a finite size (when b < 0) or infinite size (when b > 0) 
and then turns around to repeat the cycle. In these solutions, as described in detail in [l(|[ll[, ^ ne behavior 
of a% (t) near the singularity at t ~ is smooth (if po or K take generic values, then a? E (r) ~ t 2 — > 0; if 
both po,K vanish or take some special values, then a% (t) ~ t 6 — > 0). Also, the behavior of the potential 
and kinetic energy terms for the scalar field a (r) near the singularity are surprising. Namely, contrary to the 
generic solution, the potential energy dominates over the kinetic energy so that the equation of state w (r) is 
negative near the singularity. According to common lore, it was thought that such zero-size-bounce solutions 
would not exist because they would violate the null energy condition (NEC). However, this is not the case. 
The NEC is satisfied because there is a singularity at zero size, and this is the exception allowed according to 
the NEC theorems. We found all such solutions and classified them in [l(|[llj], ^ nus providing a rich class of 
examples of zero-size-bounce cyclic universes, characterized by arbitrary values of the parameters po, K, b, c plus 
one additional quantized parameter. 

To make the zero-size-bounce happen, a synchronization of initial conditions and a quantization condition 
among the 6 available parameters must be imposed. These properties are illustrated in Fig. 3, where it is 
seen s 7 (0) = (f>~ (0) = is imposed as an initial condition, and the periods of oscillations of <fi (r) and s (r) 
are quantized relative to each other. Hence such solutions are characterized by 4 continuous and 1 quantized 
parameter, rather than the 6 continuous parameters of the generic solutions. In that sense the zero-size-bounce 
solutions are a set of measure zero. So, statistically, it does not seem likely that the universe would choose such 
a solution over the generic solution. 

Fig. 4 - Geodesically complete trajectory only in Fig. 5 - Geodesically complete trajectory only in the 
the gravity sector, for b < 0. gravity sector, for b > 0. 

However, anisotropy plays a very important role in providing an attractor mechanism such that, for typical 
initial conditions away from the singularity, all trajectories are attracted to the origin of the ((f), s) plane and 
can cross the lightcone in Fig.l only at the origin. In that sense the type of solution depicted in Figs. (2-5) is 
not too far from being generic once anisotropy is taken into account. In particular the behavior away from the 
origin is a good approximation. Nevertheless, what goes on in the singularity region is quite different as seen 
in Figs. (8 ,9) and discussed below. 
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Fig. 6 - Finite-size-bounce, Fig. 7 - Finite-size-bounce, 

l^min 

(r)| > K/Ab > | 

-max 

(r)| . a% (t) never vanishes. 



The finite-size-bounce (see Figs. (6,7)) describes a universe that contracts up to a minimum non-zero size 
and then bounces back into an expansion phase up to infinite size. As the universe turns around to repeat 
such cycles the minimum size is not necessarily the same, as this depends on the parameters. Such solutions 
occur when the parameters satisfy, po < K 2 /16b and ^ in (t ) > K/4b > s^ ax (to). Note that there are still 6 
parameters, so this is not a set of measure zero, but it is a restricted region of parameter space or initial values. 
The analytic solution is given explicitly in [H[ • 



IV. ANISOTROPY AND THE ATTRACTOR 



It is useful to write the equations of motion for the fields a, a±, oti in terms of the canonical momenta instead 
of their velocities. These are given by za — p„ , za\ — p\ and zct2 = P2- Then the energy constraint that follows 



from the action S^ff or ^cs takes the following form when written in terms of z (r) = 0^ (r) 



(r) 



4z 3 



2z 



6Kv (o;i, a%) 



(15) 



where p\ = p\ + p\. For the generic solution (i.e. unlike the zero-size-bounce solution, assuming typically 



(po-, pi, pi) are non- vanishing at the singularity), the kinetic term 



2z 3 



is the dominant term in the energy 



equation as we approach the singularity z — > 0. The next to the leading term is while the potential terms 

V (a) , Kv ^' a ^ are subdominant. Keeping the leading and next to the leading terms, it is possible to integrate 
all the Friedmann equations to obtain the behavior of the general solution near the singularity. In the absence of 
the potentials, the momenta {p<j,Pi,P2) are all conserved; then the solution is uniquely given as follows (where 
we use units with k = y6) 



z(t) = It VpI+pI+Pqt 



4(r) 



|z(r)|, a(r) 



In 



{r/Tf 



(VpT 



pi 



POT 



(16) 



as well as 



h (r) + s 1 (r) = Vf (VpI+pI + Par) 



(r/TY 



(17) 



Pot 



(t) - s 7 (r) 



2r 

7t 



(r/TY 



pi 



Pot 



(18) 
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where T is an integration constant that amounts to an initial value for a. 
near the singularity are 



Similarly, the solutions for ol\ (t) , (t) 



OL{ (r) 



pi 



In 



(r/T t Y 



1,2. 



Pa 



POT 



(19) 



where Ti are integration constants that amount to initial values for the Ui. 
The solutions for (z (r) , a (r)) are plotted in Fig. 8, where we see that z (r) 



(t) — s 2 (t) changes sign. This 



shows that the trajectory of the universe inevitably passes through antigravity. Furthermore, the parametric 
plot in Fig. 9 shows that there is an attractor mechanism that forces the trajectory to pass through the origin 



of the 



) plane. 




Fig. 8 - Plots of z (r) and a (r) , with all Fig. 9 - Antigravity loop in the parametric plot of 

parameters po,p a ,p a , T, set to 1 numerically. </> (t) , s (r) . 



From these plots it is seen that the universe contracts to zero size with a big crunch, passing through the 
singularity at z (— r c ) = 0, it then expands in a region of antigravity while z (r) < 0, then contracts again to 
zero size at z (0) = 0, and then emerges with a big bang in a region of gravity where z (t) > 0. The duration of 
the antigravity phase is 

r c = M±A. (20) 
Po 

The essential parameters that control the properties of this unique solution in the vicinity of the singularity are 
Per, \Jp\ + P 2 -> Po- They are related to the kinetic energy of the sigma field, anisotropy and "radiation" (i.e. all 
relativistic particles in the early universe) respectively. They are all typically non-zero and not negligible in the 
vicinity of the singularity. Note that even if anisotropy \/p\ + p\ is infinitesimal, the attraction to the origin 
(j> = s = is inevitable. 

The parameter space for which the potentials can be neglected is estimated by inserting the solution back 
into the the equations and comparing the magnitude of the potential energy terms to the kinetic energy and 
radiation terms. From this we determined that the radiation term po/ z 2 (r) in Eq. (|15[) beats the potential V (a) 
by more than two powers of r as r — ¥ 0, and similarly for r —¥ — r c , so V (o~ (r)) is negligible near the singularity. 
Similarly Kv (oti, 02) / z is negligible compared to radiation as a function of time near the singularity (even when 
K 7 ^0) provided p\ > 15 (pf + p?,) . This condition is also compatible with a weaker condition (p 2 > 3 (pf+ g|) , 
see [Tl|) that avoids the mixmaster behavior discussed in [13| and is consistent with the discussions in |18|[12|. 

The attractor mechanism is an almost model independent universal behavior since the potential energy V (cr) 
is negligible close to the singularity. Anisotropy in even tiny amounts forces the universe to behave as described 
above. The issue of initial conditions for all degrees of freedom in our model becomes partially resolved because 
at the singularity z = we must have the following automatically enforced initial values, independent of any 
other arbitrarily chosen initial conditions at some other time, 

(f> 1 = s 1 = 0, (21) 

rather than only cf> 2 — s 2 , = z = 0. Hence, we must take </> 7 (0) = s 7 (0) = as the dynamically enforced 
boundary values at all cosmological singularities. This reduces the number of arbitrary initial values. 
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Now we see that the zero-size-bounce solutions described in the previous section (Figs. (2-5)) are not really of 
measure zero. They do take into account effectively the dynamically enforced initial values of Eq. ([2~lj) . Although 
these solutions are not capable of describing the detailed (antigravity loop) behavior near the singularity, they do 
describe the correct analytic behavior away from the singularity where anisotropy becomes negligibly small. A 
more accurate geodesically complete solution is the combination of the near-singularity and far-from-singularity 
solutions which can be constructed by matching them at some time sufficiently away from the singularity. 

V. ANTIGRAVITY REGIME 

The physics during the antigravity period remains cloudy in our investigation so far. We discussed the 
trajectory of the universe on the average, including the antigravity period. However, once in this regime, we 
may expect some violent behavior and backreaction because of antigravity, and these may affect the trajectory 
during the antigravity period. First, because of the opposite sign of the gravitational constant we expect that 
chunks of similar matter will repel each other rather than attract. This is an invitation to investigate the nature 
of interactions, including fluctuations on top of our background solution, which so far we have not included in 
the analysis. Furthermore, there are additional phenomena we should expect due to interactions as seen by 
investigating the form of the action in the antigravity regime, as follows. 

Because of the sign change of the effective Newton constant one may argue that fluctuations of the metric, 
namely the spin-2 gravitons, will now come with the wrong sign kinetic terms. There are no gauge symmetries to 
remove them, hence, they will behave like negative energy fluctuations, which could destabilize the vacuum by 
the emission of negative energy gravitons together with lots of positive energy matter particles. There is however 
a counter effect. Once particle production begins, it automatically implies a back reaction that changes the 
energy density p for all relativistic matter that enters in the cosmological equations. This effectively amounts 
to a huge increase of the parameter po once we enter the antigravity regime. So, rather than solving the 
cosmological equations with a constant po m the antigravity region, we should allow po to increase dramatically 
by large amounts. Then according to our rough formula in Eg. (1201) for the duration of the antigravity regime, the 
duration gets shortened as po increases. In other words, particle production by antigravity shortens its duration. 
The particle production and the corresponding jump of po need to be better understood. These are effects that 
oppose each other. It is difficult to estimate the balance of these effects with our current understanding of the 
corresponding physics. 

Of course, quantum gravity effects need to be also included. Therefore, it would be very interesting to 
study similar circumstances in the framework of string theory, although the requisite technical tools may not 
be available at this time. To formulate the antigravity aspects in string theory we could use the BCST field 
transformations given in Eqs. (|9lll[) . but even better would be the inclusion of the analog of the Weyl symmetry 
in the framework of string theory. 

We should also mention that the phenomenon of the sign change of the Newton constant transcends the specific 
simple model in Eg. (1131) . The phenomena we have found should also be expected generically in supergravity 
theories coupled to matter whose formulation include a similar factor that multiplies R (g) . In supergravity, 
that factor is related to what is called the Kahler potential, and this factor, combined with the usual Einstein- 
Hilbcrt term, is not generally positive definite [15| . In the past, typically it was assumed that the overall factor 
is positive and investigations of supergravity proceeded only in the positive regime. A discussion of the field 
space in the positive sector for general Af=2 supergravity can be found in [lq . However, our results suggest 
that generically the overall factor can and will change sign dynamically, and therefore antigravity sectors similar 
to our discussion in this paper should be expected in typical supergravity theories. Furthermore, we can also 
draw an analogy to similar factors that occur in front of the kinetic terms for gauge bosons and scalar fields 
[lj|. These too can potentially change sign. In the case of gauge bosons, since the kinetic term takes the 
form (E 2 — _B 2 ) where E,B represent electric and magnetic fields, a sign change could be reinterpreted by an 
electric-magnetic duality transformation. Indeed this type of duality was explored in depth at the quantum 
level for Af—2 super Yang-Mills theory in [l7|. It is therefore interesting to point out that something similar 
may be part of the phenomena in the case of the gravity/ antigravity transitions we have uncovered. 

Having mentioned a few of the possible physical effects in the antigravity regime, we will proceed by assuming 
that the destabilizing effects of antigravity are not a disaster, but rather they are the signals of interesting 
physics that remains to be understood better. We expect that the general nature of our solution survives 
and still describes the interesting effects we discovered, namely the attractor mechanism and the presence of 
antigravity for some period of time between a crunch and the following bang. 
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VI. SUMMARY 

In this talk I emphasized the role of analytic solutions of cosmological equations to discuss geodesic complete- 
ness through the big bang singularity and the discovery of some new phenomena by focussing on questions that 
could not be answered with only approximate solutions. Until better understood in the context of quantum 
gravity, or string theory, our results should be considered to be a first pass for the types of new questions they 
raise and the answers they provide. Here is a brief summary of the main points we have learned. 

• We found new techniques to solve cosmological equations analytically. This led to all the solutions for 
several special potentials V(a) (only one of those was discussed in this talk). Several model independent 
general results that followed from this include geodesic completeness, and an attractor mechanism to the 
origin of the (</> 7 , s 7 ) plane, which resolves partially questions on initial values. 

• Antigravity is very hard to avoid. Anisotropy + radiation + kinetic energy of a, require the antigravity 
epoch. For close to a year we tried to find models and mechanisms to avoid antigravity (i.e. when all 
solutions of a geodesically complete model are included). The failure to find such mechanisms finally led 
us to take antigravity seriously. 

• We studied the Wheeler-deWitt equation to take into account some quantum effects for the same system 
that we analyzed classically. We could solve some cases exactly, others semi-classically. This will be 
published at a later time. But let me mention that we arrived substantially to the same conclusions that 
we learned by studying the purely classical system. 

• Will the new insights we obtained survive the effects of a full quantum theory? My expectation is 
affirmative. Perhaps the most satisfactory approach is to try to study the same issues in the context 
of string theory. However, this is actually quite challenging since string theory has not been developed 
sufficiently in the cosmological context. 

• I emphasize that the structure of the theory in Eq .(fT"3|). the methods of solutions, and the phenomena that 
ensued, are historically direct predictions of 2T-physics in 4+2 dimensions [HQ. The model corresponds 
to the conformal shadow of 2T-physics which takes the form of a familiar lT-ficld theory setting. Other 
shadows [3| of the same 2T theory which would be in the form of dual lT-field theories, could also play 
an important role to understand the system and its physical interpretation. 

• Open questions include: are there observational effects today of a past antigravity period? This is an 
important ongoing project for our group. This involves the study of small fluctuations and the fitting to 
current observations of the CMB. 

Much remains to be understood, including quantum gravity and string theory effects, but it is clear that 
previously unsuspected phenomena, including an attractor mechanism and antigravity, are at work close to the 
cosmological singularity. The technical tools to study such issues in the context of a full quantum theory of 
gravity are yet to be developed. So this is an important challenge to the theory community, including string 
theorists, since here are some deep physics issues where string theory, if it is the correct approach at those scales, 
should play its most important role as a quantum theory of gravity. Perhaps duality transformations of the 
type that occur in electric-magnetic duality [l7j could play a role as mentioned earlier. Our understanding of 
cosmology and of the beginning of the universe could not be clarified until an understanding of the antigravity 
effects are obtained in a full quantum theory of gravity, and the associated phenomena are taken into account 
in comparing to cosmological observations. 
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